arXiv: 1508.06278v2 [hep-th] 25 Sep 2015 


Quantum Electrodynamics in d = 3 from the e-expansion 


Lorenzo Di Pietro,* Zohar Komargodski/ Itamar Shamir,* and Emmanuel Stamou§ 
Department of Particle Physics and Astrophysics 
Weizmann Institute of Science, Rehovot 7610001, Israel 
(Dated: August 26, 2015) 

We study Quantum Electrodynamics in d = 3 (QED 3 ) coupled to Nf flavors of fermions. The 
theory flows to an IR fixed point for Nf larger than some critical number Nf. For Nf < Nf, chiral- 
symmetry breaking is believed to take place. In analogy with the Wilson-Fisher description of the 
critical O(N) models in d = 3, we make use of the existence of a fixed point in d = 4 — 2e to study 
the three-dimensional conformal theory. We compute in perturbation theory the IR dimensions 
of fermion bilinear and quadrilinear operators. For small Nf, a quadrilinear operator can become 
relevant in the IR and destabilize the fixed point. Therefore, the e-expansion can be used to estimate 
Nf. An interesting novelty compared to the O(N) models is that the theory in d = 3 has an enhanced 
symmetry due to the structure of 3d spinors. We identify the operators in d = 4 — 2e that correspond 
to the additional conserved currents at d = 3 and compute their infrared dimensions. 
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I. INTRODUCTION 

We consider an K. gauge theory in d = 2 + 1 dimen¬ 
sions, coupled to 2 Nf complex two-component massless 
fermions of unit charge, (i = 1,..., 2 Nf). This theory 
has an SU(2Nf) global symmetry . 52 When Nf is suffi¬ 
ciently large, the theory flows to a stable interacting fixed 
point with an SU(2Nf) global symmetry. (It is stable in 
the sense that there are no relevant operators preserving 
all the symmetries.) However, the IR behavior is differ¬ 
ent if the number of fermions is smaller than a critical 
value Nf < Nf, leading to spontaneous symmetry break¬ 
ing according to the pattern 1 

SU{2N f ) SU(Nf) x SU{N f ) x 17(1) . (1) 

This symmetry breaking pattern can be triggered by the 
condensation of the parity-even operator 

N f 

i:^r-t> a+ N f r +N ‘) ■ ( 2 ) 

a —1 

Various estimates of the critical number Nf exist in the 
literature. 

In condensed matter physics this theory has been advo¬ 
cated as an effective description of various strongly corre¬ 
lated materials. QED 3 can arise as the continuum limit of 
spin systems with various values of Nf, e.g. Nf = 2,4 . 2-4 
The theory with Nf = 2 has also applications in high- 
temperature superconductivity . 5-7 

A method that has been employed to study QED 3 is 
the large- Nf expansion . 8-11 At large Nf, the theory sim¬ 
plifies and a systematic expansion in 1 /Nf can be carried 
out . 53 For an alternative to large Nf that uses the func¬ 
tional renormalization group approach see 14 . 

Here, we study QED 3 using the e-expansion. Clearly, 
since the theory is IR free in d = 4 and since the gauge 
coupling has positive mass dimension for d < 4, there is 
an IR fixed point at d = 4 — 2e with e > 0. The fixed 
point is generated analogously to the Wilson-Fisher fixed 
point . 15 Experience with the Wilson-Fisher fixed point 


of O(N) models, on which both e and large-IV expan¬ 
sions have been applied , 54 suggests that in QED the e- 
expansion may be more effective in describing small- Nf 
physics. 

In the development of the e-expansion for QED one en¬ 
counters some new technical difficulties that do not arise 
for O(N) models. Perhaps one reason that (to our knowl¬ 
edge) it has not been considered before is that spinor rep¬ 
resentations of the Poincare group do not behave very 
simply as a function of the number of dimensions (un¬ 
like tensor representations ). 55 Hence, it may not be ob¬ 
vious how to analytically continue to d. However, there 
appears to be no fundamental obstruction to studying 
QED with d < 4. We will keep the spinor structure that 
exists in d = 4 also in lower dimension. In lower inte¬ 
ger dimension, the representation is reducible and can be 
interpreted in terms of the existing spinor structures in 
d = 3 and d = 2. 

The fact that spinor representations are smaller in d = 

3 than in d = 4 enhances the symmetry of the theory. 
The theory in d = 3 enjoys an SU (2 Nf) global symmetry, 
while the theory in d = 4 (around which we expand) only 
an SU(Nf) x SU(Nf) symmetry. We find that in d — 

4 — 2e certain antisymmetric tensor operators, bilinear in 
the fermions, are naturally interpreted as continuations 
of the enhanced currents of the three-dimensional theory. 
This suggests that the e-expansion provides the necessary 
elements to correctly describe the theory in d = 3. 

Here, we only perform leading-order computations in 
the e-expansion of QED. Going to higher orders will be 
necessary to acquire more confidence about the accu¬ 
racy of the method, and to estimate the uncertainties . 19 
We consider bilinear and quadrilinear operators in the 
fermions. We shall see that a certain quadrilinear op¬ 
erator invariant under SU{2Nf) and parity can become 
relevant in the IR for low values of Nf, and may desta¬ 
bilize the fixed point. At leading order in e, evaluating 
the dimension naively at e = 1/2 without any resumma- 
tion leads to Nf = 2. This is consistent both with the 
F-theorem 20 and with lattice data . 21-24 (A different esti¬ 
mate that uses input from the d = 2 flavored Schwinger 
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model gives N'j = 4. 56 ) We also estimate the dimensions 
of these bilinear and quadrilinear operators at the fixed 
point for Nf > N c f . 

A. Compact vs Non-Compact Gauge Group 

In this note we do not study QED 3 with compact 
gauge group 17(1), 25 but we would like to make sev¬ 
eral comments nevertheless. 5 ' The compact theory has 
monopole operators that sit in various representations 
of SU(2Nf) x U(1)t, where U(1)t is associated to the 
topologically conserved current j lt = ep Up F up . Monopole 
operators can condense (i.e. have an expectation value 
in the vacuum) and they can proliferate (meaning that 
we can add them to the action). Because monopole op¬ 
erators are charged under the global symmetry group, 
the question of monopole proliferation is sensitive to the 
symmetries preserved by the regulator (be it the lattice or 
some other UV theory that flows to QED 3 ). Therefore, 
in QED 3 monopole proliferation has presumably noth¬ 
ing to do with chiral-symmetry breaking, which is an 
intrinsic phenomenon, given that it also happens in the 
non-compact theory, which has no monopoles. In fact, 
if monopole operators proliferate, the Nambu-Goldstone 
bosons of ( 1 ) would generically become massive because 
the SU(2Nf) symmetry would be broken explicitly. This 
mechanism of gapping the theory is analogous to the 
mechanism for confinement in Polyakov’s model . 26 On 
the other hand, we do expect monopole operators to con¬ 
dense. The lightest monopole operator sits in the rank- 
Nf antisymmetric representation, and its condensation 
would be precisely consistent with the symmetry break¬ 
ing pattern (l ). 58 

II. THE e-EXPANSION 
A. Generalities 

To illustrate the procedure of the e-expansion, con¬ 
sider the two-point function of an operator O in d = 4, 
expanded in perturbation theory in a classically marginal 
coupling g 

(0(p)0(-p)) = p 2A ~ 4 ^2 Cnm g U ( log 2 ) > (3) 

0<m<n,n—0 ' P J 

where A is the dimension of O in d = 4 at g = 0, and 
A is a UV cutoff. Introducing the renormalized operator 
O ren = Z0 1 we can cancel the A-dependence of the cor¬ 
relator by allowing the coupling g and the normalization 
z t0 evolve according to = P{g), = 7 (g), 

such that the Callan-Symanzik equation holds 

(al^A +P{9) §g + 27(5) ) W)°(-P)> = 0 • ( 4 ) 

The terms in (3) with coefficients c nn , n > 1, are the 
leading logs. It follows from (4) that they are all fixed in 


terms of the coefficients pi and 71 in the leading order 
expansion of /3 and 7 

P{g) = Pig 2 + 0{g 3 ), 7 (g) = 7 id + 0(g 2 ) . (5) 

One can then resum the leading logs to obtain 

_ 2 ’7'1 

{0(p)0(-p)) ~ p 2A “ 4 ^1 + |/3i g log ^ . ( 6 ) 

For d = 4 — 2e we assume that g acquires a positive 
mass dimension ce (where c is some positive number). 
The analogue perturbative expansion of the two-point 
function in d = 4 — 2e is 

ww-pii-A'Ei. |-.) ■ < ? ) 

n—0 7 

Requiring that (7) approaches (3) in the limit e —» 0, we 
find the matching condition 


Cn 


n 

T. Cnm m\ 
m —0 



+ 0(e) . 


( 8 ) 


The leading contribution to the two-point function (7) 
in the limit eCl comes from the terms containing c nn , 
which we can thus resum similarly to ( 6 ) 

_ 2 T 1 

(o(p)o(~p)) ^ P ^- d ( i+^-y p 2A -v 71 ^. 

\ ce p ce J p->o 

(9) 

In the IR limit p —> 0 a new scaling law emerges. The 
contribution to the IR dimension of the operator at first 
order in e is thus 

Air = A + 71 —+ 0(e 2 ). (10) 

Pi 

The crossover to the IR scaling in (9) happens when 


1 < 


Pi g_ 

ce p ce 



( 11 ) 


We see here the physical consequence of introducing the 
small parameter e: the crossover towards the IR happens 
at a scale that is enhanced by the parametrically large 

_ 1 _ 1 

factor (/3i/ce) ce with respect to the naive scale . As 
a result, the IR fixed point is parametrically close to the 
one in the UV. 

Indeed, the IR fixed point corresponds to a zero of the 
/3-function for the dimensionless combination g = gA~ ce 

- 0(s) = - ce 9 + Pig 2 + 0(g 3 ) 

P{g*) = 0 => ()* = — + 0 (e 2 ) . ( 12 ) 

Pi 

Comparison with (10) shows explicitly that, at leading 
order in e, the difference Air — A is the anomalous di¬ 
mension 7 evaluated at the fixed point . 2 ' -28 Extrapolat¬ 
ing the results to e = 1 , we obtain an estimate for the 
observables of the IR theory in three dimensions. 
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B. Wilson-Fisher Fixed Point in QED 

The Lagrangian for QED in d = 4 is 
1 Nf 

c = + i • (13) 

a —1 

We use the usual four-dimensional Dirac notation for the 
spinors. Their decomposition in terms of two-component 
fermions is 

*° = ( ia 2 ip a+N f ) ’ Q = '' ’ ’ Nf ' ( ' 14 ^ 


Nf > Np the full SU(2Nf) is realized linearly at the 
IR fixed point. This entails the existence of 2N 2 + 1 
additional conserved operators of spin 1 with protected 
dimension A = 2. It is natural to ask whether these 
operators are visible also in the theory continued to non¬ 
integer dimension. 

One of the additional currents is the singlet axial cur¬ 
rent J® = Indeed, the continuation of the 

anomaly operator F A F vanishes for d = 3. As for the 
remaining 2Nj currents, we note that in d = 4 — 2e we 
can define the following antisymmetric tensor operators 

(K„ v ) a b = . *a7^75* b ■ (17) 


In dimension d we take the Clifford algebra to be 
= 2r]^ u t, with rj^r)^ = d. To leading non¬ 
trivial order, the beta function in d = 4 — 2e is given by 
(here e = eA” e ) 

AT, 

/3(e) = -ee + j^e 3 + 0(e 5 ) . (15) 


The value of the coupling at the Wilson-Fisher fixed point 
is e 2 = 127 r 2 e/Nf. The theory is therefore weakly coupled 
when we are close to d = 4 or when the number of flavors 
is large. 

A comment on 75 is in order. A consistent defini¬ 
tion of 75 in non-integer dimension is due to ’t Hooft 
and Veltman . 29-31 According to this prescription, 75 an¬ 
ticommutes only with the y^’s of the four-dimensional 
subspace, and commutes with all others. This implies 
an explicit breaking of axial symmetries in d = 4 — 2e, 
and reproduces the chiral anomaly for the singlet axial 
current X) a (f , a 7 Al 75 v I' a in the limit e —> 0 . 32,33 For the 
leading-order calculations that we present here, this pre¬ 
scription is in practice equivalent to a naive continuation 
of 75 as totally anticommuting. However, the difference 
from the naive continuation becomes relevant at higher 
orders. 

QED in d = 4 has an SU(Nf) x SU(Nf) global sym¬ 
metry with associated conserved currents 


Wa = *a7„* 6 - J2 ^7„* C , 

N f c 

(4)a = *a7„75’J'" - ~^S a b ^75-T • 

7 r 


(16) 


Their anomalous dimension at one-loop vanishes and, 
therefore, at leading order the IR dimension is the same 
as the classical one, i.e. d—1. This is the correct scaling 
dimension for conserved currents. For the vector current 
this argument is valid at all orders in perturbation the¬ 
ory, because they are conserved for any d. On the other 
hand, the axial currents J 3 are explicitly broken for non¬ 
integer d 32 and this can affect the IR dimension at higher 
orders. Nevertheless, we do expect them to be conserved 
in d = 3, because the non-conservation is given by an 
operator that vanishes both in d = 4 and d = 3. 

So far, we have argued that the e-expansion predicts 
the existence of currents associated to the global symme¬ 
try SU(Nf) x SU(Nf) in the IR CFT for d = 3. How¬ 
ever, QED 3 has an enhanced SU{2Nf) symmetry. For 


They carry the correct flavor and Lorentz quantum num¬ 
bers to be identified with the additional currents, because 
in d = 3 we can use the totally antisymmetric tensor e M „p 
and dualize them to spin 1 operators. 

We are led to the expectation that the IR dimension 
of J® and should evaluate to 2 for e = 1/2. The 
one-loop computation gives 

A m (J s ) =3-2e + 0(e 2 ), 

Air(AT) = 3 - 2e + + 0(e 2 ). 

The anomalous dimension of J s only starts at two- 
loop order. 33,34 As we will show in the next section, we 
can estimate that the IR critical point exists only for 
Nf > 3. Plugging Nf = 3 and e = 1/2 in (18) we find 
Aj^ oop (AQ = 2.25, which agrees with the expectation 
within a 10% margin. The precision improves for larger 
values of Nf . We view this as a hint that the continuation 
to non-integer dimensions correctly captures the proper¬ 
ties of the 3d CFT that we ultimately want to study. A 
preliminary check of higher orders in e shows that the 
agreement improves. This will be discussed elsewhere. 

Note that we can also study the anomalous dimension 
of the operator F IJV . The Bianchi identity is obeyed for 
all d, and one can verify that Air(F) = 2 holds to all 
orders in e. 


C. Quadrilinear and Bilinear Operators 

In the three-dimensional theory there are two parity- 
even quadrilinear scalar operators that are invariant un¬ 
der the full SU(2N f ) 

0 i = (^^cdV ) 2 and O 2 = (X ^^ ? ) 2 • ( 19 ) 

i i 

The operator 0\ can be easily continued to d = 4 — 2e. In 
Dirac notation we can rewrite it as 0\ = (]T a ^ a l^ a ) 2 - 
To continue the operator 0 2 , we use the fact that in 
d = 3 the antisymmetrization of three 7 matrices is pro¬ 
portional to the identity. Therefore, we can rewrite it 
as 0 2 = 6 (]C a 4'o7[^7^7p] v k a ) 2 , which is a well-defined 
expression also for d = 4 — 2e. Note that in d = 4 this 
operator can be identified with the square of the axial 
current (X) n 
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To obtain their IR dimension we compute the mix¬ 
ing between these two operators. Typical diagrams at 
one-loop are shown in Fig 1. To obtain the correct 



FIG. 1: Diagrams giving the mixing matrix of the quadrilinear 
operators at one-loop. 


mixing at one-loop, it is necessary to also take into ac¬ 
count the one-loop mixing with the operator Oeom = 
(E a *a'V li * a )G&'F liV - E fe ^ 7 ^ & ) that vanishes on 
the equations of motion. This mixing is induced by the 
diagrams in Fig 2. 



FIG. 2: Diagrams giving the mixing matrix of 0\ and O 2 into 
Geom- 


In the basis {0 1 ,0 2 }, the matrix of anomalous dimen¬ 
sions reads 35 




( 20 ) 


Its eigenvalues are j §^ j (2 Nf + 1 ± 2 y/fV 2 + Nf + 25). 59 

By evaluating them for the fixed-point value e 2 = 
127 r 2 e/Nf, we find that the operator corresponding to 
the negative eigenvalue becomes relevant when Nf < 
Y + 0(e 2 ). Thus, for e = \ the operator is relevant 
in the IR when Nf = 1,2, while it remains irrelevant for 
all integer Nf > 2. Because the quadrilinear is neutral 
under all the global symmetries, when it becomes rele¬ 
vant it may be generated and trigger a flow to a new IR 
phase (e.g. a Goldstone phase). From this we obtain the 
estimate Nj < 2. 

The same mechanism for the onset of chiral- 
synnnetry breaking has been studied using the large- Nf 
expansion . 36 The large- Nf approximation of the anoma¬ 
lous dimensions is such that 0 1 and O 2 are always irrel¬ 
evant at the IR fixed point for Nf > 1. (See however the 
renormalization group study at large Nf in 37 .) 

Let us make a few comments on d — 2. There, the 
quadrilinear operator is marginal already at the tree- 
level. Therefore, the criterion above implies that in d = 2 
the anomalous dimension must evaluate to 0 at Nf. This 
is satisfied for Nj = 00 . This value is consistent with the 


IR behavior of QED in d = 2: for every finite Nf the the¬ 
ory flows to an SU(2Nf) Wess-Zumino-Witten (WZW) 
interacting CFT , 38,39 namely a cr-model with coset target 
space deformed by a WZW term (exactly as one expects 
for the Lagrangian of Nanrbu-Goldstone bosons). Even 
though Nambu-Goldstone bosons do not exist in d = 2, it 
appears natural to interpret this theory as the continua¬ 
tion of the cliirally broken phase to d = 2. See also 40,41 . 
Assuming that the anomalous dimension approaches 0 as 
1/Nf for Nf —► 00 , the divergence of Nj(d) for d = 2 is 
given by a simple pole. This suggests using the modi¬ 
fied ansatz Nj{d) = (d — 2)~ 1 f(d) in the equation for 
Nj, which can then be solved for / perturbatively in 
e. With this ansatz, the leading-order estimate becomes 
Nj < 4.5. The difference with the previous estimate are 
liigher-order terms in e; it may be viewed as a measure 
of uncertainty. Improving on this requires computations 
beyond one loop . 19 

Further data about the fixed point can be obtained 
by considering bilinear scalar operators. There are two 
types of scalar operators in the three-dimensional theory. 
Operators of the first type are scalars also in d ^ 3, i.e. 

(5i) 0 b = T q 4A , $ a7 5 ^ . (21) 

Operators in this class preserve at most the diagonal 
SU(Nf ) subgroup of SU(Nf ) x SU(Nf), and the most 
symmetric ones are ^2 a ip a ' l l ,a+Nf A c.c.. The one-loop 
computation (see Fig 3) gives 

A m {B 1 )=3-2e- — + G(e 2 ). ( 22 ) 



FIG. 3: The diagram giving the anomalous dimension of bi¬ 
linear operators at one-loop. 


The second type of scalar operators are given by rank- 
three antisymmetric tensors in d = 4 — 2e 

{B2^ P )a = ’f'a7[ M 7.7p]^ b , ^a7[p7.7p]7 5 ^ 6 • (23) 

They give rise to scalars in d = 3 because they can be 
contracted with the totally antisymmetric tensor e M1/p . 
The chiral condensate (2) and the parity-odd, SU{2Nf)~ 
invariant bilinear belong to this class of opera¬ 

tors. Since their anomalous dimension vanishes at lead¬ 
ing order in perturbation theory, their IR dimension to 
first order in e is captured by just the classical contribu¬ 
tion 

Air(B 2 ) = 3 — 2 e + 0 (e 2 ) . (24) 

The anomalous dimension starts being non-zero at two- 
loop order . 33,42 This implies that higher orders in e in 
(24) will be non-zero. Nevertheless, the e-expansion sug¬ 
gests that the IR dimension of these scalar operators is 
perhaps close to A = 2. 
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III. FUTURE DIRECTIONS 

In this paper we initiated a study of the critical point of 
QED 3 based on the e-expansion. Our results were based 
on leading-order computations. It would be very inter¬ 
esting to sharpen the theoretical predictions by higher- 
order computations . 60 The necessary preliminary step 
of computing the two-loops counterterms was done for 
generic gauge theories with fermions in 45 . Due to the 
asymptotic nature of the e-expansion, efficiently includ¬ 
ing higher-order terms requires the use of resummation 
techniques. (In this context, it would be interesting to 
understand if data from the flavored Schwinger model 
can be efficiently included.) 

The e-expansion can be used to compute additional ob¬ 
servables of the IR fixed point. For instance, correlators 
of the stress tensor and of conserved currents . 61 Another 
interesting datum of the 3d theory is the universal co¬ 
efficient F of the partition function on the three-sphere, 
which gives also the universal part of the entanglement 
entropy across a circular region. For this one can utilize 
the techniques of 48,49 . The calculation of F in QED 3 


via the e-expansion has been recently presented in 50 . 
Another interesting line of investigation would be to see 
if the conformal bootstrap techniques shed any light on 
QED 3 . 51 
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